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ABSTRACT 


This  study  was  undertaken  in  order  to  extend 
previously  developed  computational  methods  for 
conventional  dome-transducer  interactions  to  the 
nonlinear  pare  metric  case.  A  simple  model  involving 
planar  geometry  was  analyzed  in  detail  in  order  to 
gain  insight  into  the  sensitivity  of  parametric 
generation  efficiency  to  variations  in  such  param¬ 
eters  as  steering  angle,  dome -transducer  spacing, 
and  dome  thickness.  A  prescription  is  given  for 
calculating  the  interaction  for  realistic  geometries. 
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INTRODUCTION 


1.0 

In  the  past  Tracor'*'  has  developed  computer  programs 
which  allow  predictions  to  be  made  of  the  influence  of  a  dome 
upon  the  radiated  field  patterns  of  a  conventional  sonar. 
Programs  were  developed  to  cope  with  a  variety  of  realistic 
geometrical  configurations.  Moreover,  formulations  were  given 
for  thick  shells  which  included  structure.  All  of  these  were 
based  upon  small  signal  theory. 

Thu  objective  of  the  current  program  has  been  to 
extend  our  computational  ability  for  the  dome- transducer  inter¬ 
action  to  the  nonlinear  parametric  case. 

Simplifying  assumptions  are  given  along  with  their 
justification.  We  have  utilized  these  assumptions  to  formulate 
a  model  in  the  idealistically  simple  case  of  an  infinite  planar 
dome -transducer  interaction.  From  this  model  we  can  learn 
something  about  the  sensitivity  of  parametric  generation  effi¬ 
ciency  to  variations  in  steering  angle,  dome  thickness,  and 
dome  material  properties.  We  can  not  learn  anything  about 
variation  in  radiation  patterns  from  such  a  simple  model.  None¬ 
theless,  by  studying  the  planar  case  we  feel  that  we  have  gain¬ 
ed  sufficient  confidence  to  be  able  to  recommend  a  procedure 
for  computing  the  dome-transducer  interaction  for  more  real¬ 
istic  geometries  drawing  upon  previously  developed  computer 
programs . 
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la  the  next  "wo  sections,  the  derivation^  of  our 
planar  model  will  be  presented.  A  subsequent  section  will  give 
numerical  results  for  both  thick  and  thin  pl.anar  domes  showing 
the  transmission  loss. 

2.0  PRIMARY  FIELD  SOLUTION* 

The  infinite  planar  dome- transducer  model  is  shewn 
in  Fig.  1,  along  with  the  coordinate  geometry.  A  thin  plate 
is  located  a  distance  a  from  the  planar  "transducer”.  The 
transducer  surface  has  a  continuous  normal  velocity  distribu¬ 
tion  which  simulates  the  particle  velocity  of  a  plane  wave 
traveling  in  a  direction  making  an  angle  0  with  the  normal 
to  the  plane.  This  velocity  distribution  is  given  by 


(1) 


where  V0  is  just  n  velocity  amplitude  factor.  The  pressure 
field  between  the  dome  and  transducer  will  be  denoted  by  p^, 
and  that  outside  the  dome  by  p^.  In  order  for  the  fluid 
particle  velocity  to  equal  the  transducer  velocity  at  y  =  0, 


__1  op  ! 
iuj p  oyj 

y^O 


V. 


(2) 
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The  plate  is  assumed  to  be  thin  enough  so  that  the 
boundary  conditions  at  its  surface®  may  be  matched  at  the 
middle  surface.  One  condition  is  that  the  fluid  particle 
velocity  at  the  inner  and  outer  plate  surfaces  equal  the 
plate  velocity .  If  the  plate  displacement  from  equilibrium 
is  denoted  by  u,  the  equality  of  velocities  is  given  by 


iM,  lr35?2,  £u  (3) 

y„d  ^sT'-'y-d  “  snr  • 

The  dynamic  behavior  of  the  plate  is  described  by  the  classical 
plate  equation.  In  this  case, the  driving  force  is  the  pressure 
difference  across  the  plate,  so  that  the  equation  of  motion  is 


Cpx  -  p23 


a  2  ..  4 

_  u  o  u  ,  r.h  o  u 
P  cn  — *•  + - *  — 7- 

s  dtz  12(l-v2)  dx4 


Here  pg  is  the  plate  mass  density,  h  is  its  thickness,  E  is 
Young's  modulus,  and  *  is  Poisson's  ratio. 

The  solutions  to  the  wave  equation  for  the  pressure 
fields  inside  and  outside  of  the  dome  will  have  the  forms 


Pl  =  [A  eiT,y  +  B  e-iriy]eiKx  ,  r,2  +  K2  -  k2  (5) 

P2  -  C  eiKx  .  (6) 
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Substituting  these  equations  into  the  boundary  condition 
equations  yields  expressions  for  the  unknown  coefficients  in 

Eqs.  (4)  and  (5).  The  final  results  for  the  pressure  fields 
are 


Pi 


IiDp  V  .  .  ,  .. 

D  le  ^  y  +  -f  HK4)  cos  tiCy-d)}^*^* 

U)  p  - 


p,  -  iJJk  eiTr(y-d)  eiKx 


(7) 

(8) 


where  cr  »  p  h  ,  K  » 


Eh* 


12(l-v*) 

D  *  in{e“1Tld  +  -y~(-aoj2  *«-  HK4)  sin  ^d]  . 


2  ,  K  -  k  sin  6  ,  rj  «  k  cos  9 


u»  p 


It  is  seen  that  is  the  sum  of  the  traveling  wave  which  is 
transmitted  through  the  dome  and  a  standing  wave  which  always 
has  an  antinode  at  the  plate.  The  average  intensity  passing 
through  the  plate  is  just 


(upV0)‘ 


o  2oc 


W 


(9) 


For  many  cases  cf  interest  the  HK4  term  in  D 

o 

is  small  relative  to  <rw  .  and  we  may  write 


D  «  in  |e"iTld  -  h  Is. 


h  rj  sin  -ndj. 


-5- 


Defining  T  as  the  ratio  of  ^  in  Eq.  (9)  to  the  value  that  ^ 

would  have  with  no  dome  'i.e. ,  h  »  0),  we  have 

•  P  P  i 

T  »  1/jl  +  Tih  ~  sin  r,d  [nh  ~  sin  -qd  -  2  cos  nd'jj-.  (10) 

It  is  seen  that  T  is  a  periodic  function  of  d,  with  period 
▼/k  cos  9  ,  and  that  its  amplitude  of  oscil1ation  decreases 
as  6  increases.  Starting  from  d  ®  0,  T  increases  from  unity 
until  its  first,  maximum  occurs  at  a  value  of  d  for  which 

tan  2  r,d  =  2/r\h  ps/p  .  (11) 

O 

Typically,  ps/p  «7.9  and  kh  w  l.b  x  10  ,  so  that,  for  9  **  0, 

Eq.  (II)  becomes  approximately 

tan  2  *d  «  14.1  . 


This  implies  a  plate- transducer  spacing  of  .12  wavelengths. 
Maxima  occur  at  this  value  of  d  and  at  every  half  wavelength 
thereafter. 


To  extend  this  solution  to  parametric  excitation 
of  the  transducer,  the  velocity  distribution  is  given  by 


(12) 


Since  the  primary  fields  given  in  Eqs.  (7)  and  (8) 
are  obtained  by  solving  the  linear  Helmholtz  equation  then 
superposition  holds  and  we  can  immediately  write  down  the 
solutions  corresponding  to  Eq.  (12).  For  example,  in  region 
2  we  shall  have 
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ipV0 


(y-d) 

e 


iKlx 

e 


i^y-d) 

e 


iK2x 

e 


(13) 


2.0  SECONDARY  FIELD  SOLUTION 

The  secondary  radiation  field  ?_  is  given  by  the 
linearized  Helmholtz  equation 


A  a* 
?rt  P 

JJL 


(14) 


3.8  for  sea  water;  P^,  P2  are  the  time  independent  primary 
field  solutions  outside  the  dome  at  the  two  primary  frequencies. 
The  solution  of  Eq.  (14)  is  given  by 


2 

k_  P 


Poo  P-rco 


r* 


A  A* 

?2  P2  dxdy 


(15) 


In  practice,  the  00  limit  on  y  can  be  replaced  by  a  finite 

A  •  t 

limiting  distance  when  P9?2  becomes  much  less  than  unity. 
Inside  the  dome  the  solution  is 
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'gRgcORt 


A  A* 

P1  P1  * 


(16) 


I?l*  P^  are  the  time  independent  primary  field  solutions  inside 
the  dome.  Now  and  ^  are  the  Green's  functions  for  inside 
and  outside  the  dome  respectively. 

2 

An  appropriate  Gre^u's  function  for  region  1  inside 
the  dome  is  given  by  the  method  of  images  as  ' 


r>  =  -<u 


where 


r'  «  a  j2nd  +  y,  4-  a  x 
t)  y  t  J  x  o 

r  =  a  2nd  -  ya  +  a  x„ 

r,  y  L  X  o 


In  region  1  inside 


and 


the  dome,  the  boundary  conditions  on  are: 
0  at  y  =  0  (transducer) 

g  at  y  «  d 


where  g  in  most  cases  will  oe  much  less  than  unity. 

In  regicn  ?  outside  the  dome  must  also  satisfy 
the  condition  r  «  g.  For  many  applications  the  dome- transducer 
spacing  will  only  be  a  few  wavelengths  of  the  primary  beams. 

To  a  good  approximation  no  secondary  field  has  been  generated 
at  the  dome.  Therefore,  the  boimdary  condition  at  the  dome 
becomes  -  0. 

An  appropriate  Green’s  function  for  region  2  is 


(■I') 


*•1 


H, 


w/t 

1  X 


:.r)  -  H„  a)(k.R0) 


.  (18) 


Fig.  2  shows  R  and  R0  for  a  source  point  outside  the  dome. 
The  structure  of  Eqs.  (17)  and  (18)  already  guarantees  that 
r2  vanishes  at  r  »  o>  and  that  both  and  are  singular 
st  r  =  ro  . 

Without  greatly  complicating  Eqs.  (15)  and  (16), 
we  can  allow  either  for  small  signal  attenuation  or  finite 
amplitude  saturation  effects. 

The  following  parameters  can  be  varied: 

1.  primary  frequencies 

2.  dome  elastic  properties 

3.  -dome-transducer  spacing 

4.  angle  of  incidence  of  beams  on  dome 
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TRANSDUCER 


FIG.  2  GEOMETRY  PERTINENT  TO  AND 
THE  CALCULATIO'4  OF  R>  (»-). 
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If  it  is  desirable,  the  active  portion  of  the  plane  transducer 
can  be  limited  along  the  x  coordinate.  This  is  accomplished 
by  replacing  the  continuous  normal  velocity  distribution  of 
Eq.  (1)  by 


V  =  V0eikx  Sin0 
-  0 


a  . 


(19) 


The  velocity  distribution  V  given  by  Eq.  (19)  can  be  represented 
by  a  spatial  fourier  series  in  an  approximate  fashion. 

4.0  THE  INTEGRAL  FOR  P2  (  «.) 

If  we  ignore  any  generation  of  P-^(w_)  in  region  1 
then  satisfies  the  boundary  condition  ^  =  0  at  the  dome 
surface.  As  shown  in  Fig.  2  it  is  perfectly  consistent  with 
this  assumption  to  consider  our  origin  of  coordinates  on 
the  done  surface.  With  reference  to  Fig.  2,  we  see  that 


2  2  2 

R  »  r0  r  -  2r0r  cos 

(4>-</') 

(20) 

and 

2  2  2 

R0  =  ro  +  r  +  2r0  r  cos 

r—t 

CM 

v-/ 

For  the  infinite  planar  geometry  of  our  problem,  it  is  not 
correct  to  distinguish  between  near  and  far  fields.  The 
virtual  parametric  array  is  theoretically  a  tapered  end 

fire  array  that  is  infinitely  long  (  r0  - p-  oo  ).  Practically 

speaking,  it  has  significant  source  strength  for  a  finite 

distance.  For  our  purposes,  it  will  be  permissible  to 

i<  1 

truncate  this  array  at  r0  «  •  ■  -----  where  a  and  are  the 

12 
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small  signal  attenuation  coefficients  of  the  primary  beams.  If 
we  take  our  field  point  r  to  be  ten  times  larger  than  the  maui- 

k 

mum  value  for  r0  ,  i.e.,  ro  ,  then  Eqs.  (20)  and  (21)  can  be 
simplified  and  take  the  form 


R  =>  r  -  r0  cos  (<t>~4') 


(22) 


and  Ro  =  r  +  r0  cos  (4>-vf0  .  (23) 

Moreover,  we  can  choose  r  large  enough,  so  that 
k-r  >>1.  Then  the  .  arge  argument  assymptotic  formula  for  the 
Hankel  function  can  be  used  to  simplify  ^  Thus  becomes 


If  we  substitute  Eq.  (24)  into  Eq.  (15)  and  change 
to  polar  coordinates  we  obtain 


-  i  k.  r5  sin  (4+6) 

*  e  r0dr0  ^  • 

(25) 
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MMJgi/ 

Here  A  is  given  by 


3/2  f 

0Pek-  2r  * 


+  die.  cos  $  j 


A  _  * 

pc  D, 

**  0  ulu2 


(26) 


and  2  represents  the  attenuation  path  length  of  the  primary 
beams  to  the  source  point  as  shown  in  Fig.  2. 


r0  s  in  ^ 
cos0 


(27) 


Eq.  (25)  allows  us  to  calculate  the  amplitude  of  the 
difference  frequency  field  at  any  point  (r,  4> ) .  It  is  to  be 
noted  that  the  physical  parameters  relating  to  the  thin  dome 
enter  only  into  the  coefficient  A  via  and  .  In  the 
thick  dome  case,  this  coefficient  will  also  contain  the  dome 

parameters.  By  letting  the  dome  thickness  h - >  0,  we 

can  find  the  limiting  value  for  A  in  the  absence  of  a  dome. 

If  we  are  only  interested  in  the  relative  change  in  amplitude 
of  P.  at  a  field  point  (r,  b),  then  it  is  onP  necessary  to 
compute  the  following  ratio 


A  (with  dome  thick  or  thin  moc’el) 
limit  A  (thin) 

h — >0  .  (28) 

The  integral  in  Eq.  (25)  is  the  same  for  both  numerator  and 
denominator  of  Eq.  (28)  and  can  be  cancelled. 
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The  expression  for  A  in  the  thick  dome  case  was  de- 

3 

rived  with  results  given  by  Brekhovskikh  .  A  detailed  derivation 
was  presented  in  the  Third  Quarterly  Report,  and  will  not  be 
repeated  here . 

Before  presenting  numerical  results,  we  should  like 
to  list  the  key  assumptions  of  cur  planar  dome  model: 

1.  In  practice,  the  dome- transducer  spacing  will 
only  be  of  the  order  of  10  primary  wavelengths.  Therefore  all 
parametric  generation  both  within  the  dome  itself  and  between 
the  dome  and  transducer  can  be  ignored. 

2.  Small  signal  theory  is  used  to  compute  the 
transmitted  amplitudes  of  the  primary  beams. 

3.  In  a  finite  amplitude  radiation  field,  the 
finite  pressure  gradients  can  give  rise  to  hydrodynamic  flow 
termed  streaming.  Owing  to  the  proximity  of  the  dome  inner 
surface  it  is  felt  that  it  is  unlikely  that  pressure  gradients 
which  are  large  enough  to  produce  streaming  v’  11  occur.  More¬ 
over,  pulsed  source  operation  seems  to  prevent  streaming. 

4.  We  assume  that  most  of  the  energy  is  transmitted 
through  the  dome.  Eqs .  (10)  and  (11)  indicate  that  there  are 
periodic  maxima  in  the  transmission  coefficient  as  d  is  varied. 

As  will  be  seen'  in  the  next  section,  some  standing  wave  reso¬ 
nance  effects  will  occur  for  particular  values  of  d.  However 
even  for  these  cases  most  of  the  acoustical  energy  is  trans¬ 
mitted  through  the  dome;  hence,  it  is  felt  that  assumptions 

1  and  2  will  still  be  valid. 
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NUMERICAL  RESULTS 
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5.0 

In  an  effort  to  gain  some  insight  into  the  dependence 
of  parametric  array  transmitted  power  on  dome- transducer  inter¬ 
actions,  the  simple  model  of  an  infinite  planar  transducer 
and  planar  dome  was  analyzed.  While  the  infl.nite  planar  dome- 
transducer  model  cannot  be  used  to  predict  beam  patterns  and 
source  levels,  it  clearly  displays  the  dependence  of  trans¬ 
mitted  power  on  both  steering  angle  and  dome- transducer  spacing. 
Further,  the  mathematical  form  of  the  pressure  field  solution 
for  the  primary  frequencies  can  be  readily  analyzed  to  separate 
the  standing  wave  between  dome- transducer  and  the  wave  trans¬ 
mitted  through  the  dome. 

In  our  model,  the  "transducer”  was  an  infinite  plane 
with  a  continuous  normal  velocity  distribution  which  simulated 
the  particle  velocity  of  a  plane  wave  steered  in  a  direction 
making  an  angle  6  with  the  normal  to  the  plane.  The  wave 
outside  the  dome  was  just  such  a  plane  wave;  however,  the  dome 
modified  its  amplitude  and  phase.  Inside  the  dome- transducer 
cavity,  the  pressure  field  consisted  of  a  plane  wave  traveling 
in  the  outward  direction  and  a  reflected  wave  traveling  in  the 
inward  direction.  These  two  waves  included  the  direct  trans¬ 
mission  and  all  of  the  multiply-reflected  components.  The  sum 
of  all  these  components  add  to  yield  a  net  transmitted  wave 
plus  a  standing  wave. 

The  dome  was  modeled  as  either  a  lossless  elastic 

1  3 

plate  or  a  lossless  elastic  medium  .  These  two  dome  models 

are  referred  co  as  "thin”  and  "thick”,  respectively.  It  should 

be  noted  that  the  thin  model  did  not  admit  waves  internal  to  the 

dome  while  the  thick  model  did. 
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For  both  models  the  time  average  intensity  trans¬ 
mitted  through  the  plate  was  periodic  with  k  dcos  d  ,  where  k 
»  P  P 

was  one  of  the  primary  frequency  wave  numbers,  d  was  the  dome- 

transducer  spacing,  and  d  was  the  steering  angle.  This  period¬ 
icity  with  kpdcos  0  ,  coupled  with  the  fact  that  the  models  are 
lossless  and  the  plates  infinite  in  extent,  means  that  resonances 
will  exhibit  themselves  with  respect  to  both  the  dome- transducer 

•k 

cavity  and,  for  the  thick  model,  the  dome  itself.  These  effects 
can  result  in  parameter  combinations  where  the  intensity  with 
d>me,  normalized  to  the  no-dome  intensity,  exceeds  one.  A 
discussion  cf  this  for  the  linear,  single  frequency  case  can 
be  found  in  Section  2.0. 


Representative  results  of  this  analysis  are  presen¬ 
ted  in  Figs.  3  through  14.  The  parameters  used  in  this  study 
are  given  in  Table  I.  No  special  significance  should  be  attached 
to  the  parameters  selected  except  that  they  seemed  "reasonable” 
to  assess  the  effects  of  the  two  dome  models  on  the  acoustic 
interaction.  One  moderating  factor  was  that  the  only  measured 
data  we  were  aware  of  was  a  measurement  made  by  Konrad,  NUSC/NL, 
involving  a  steel  plate.  That  comparison  is  shown  in 
Fig.  15  and  will  be  discussed  in  detail  later. 


The  significance  of  resonance  is  that  we  are  changing  the 
impedance  the  transducer  sees.  Because  we  assumed  a  velocity 
controi  situation,  we  therefore  change  the  power  out  of  the 
transducer  as  a  function  of  dome-transducer  spacing  and  dome 
thickness . 

•irk 

These  results  were  communicated  privately  to  Dr.  Cary 
by  Mr.  Konrad. 
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TABLE  I 


SUMMARY,  MODEL  PARAMETERS 


Primary  Frequencies 


Dome  Material 
h 


Fluid 

c 


c 


Dome- Transducer  Spacing 
Steering  Angles 


50  kHz 
60  kHz 


Steel 


0.0  -  1.2  inches 


7.35(10) 


_4  lf>£-sec‘ 


xn. 


30(10)6  ibf/in. 2 


0.3 

2.34(10)6  in. /sec 
1.25(10)5  in. /sec 


Water 

60,000  in. /sec 

2 

-  lb --sec 
9. 66(10)"'3 

in. 

9.25  -  10.5  in. 

0°,  20° 
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NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  9.25  in.,  STEERING  0° 


NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS. 
DOME  THICKNESS,  SPACING  9.5  in.,  STEERING  0° 


MALI 2 HD  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  9.75  in.,  STEERING  0° 


NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS. 
DOME  THICKNESS,  SPACING  10.0  in.,  STEERING  0° 


TRACOR 


NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  9.5  in.,  STEERING  20° 


NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  9.75  in.,  STEERING  20° 
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NORMALIZED  FARFT.ELL)  DIFFERENCE  FREQUENCY  PRESSU 
DOME  THICKNESS,  SPACING  10.0  in.,  STEERING  2 


NORMALIZED  FARFIELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  10.25  in.,  STEERING  20° 


NORMALIZED  FAR FI ELD  DIFFERENCE  FREQUENCY  PRESSURE  VS 
DOME  THICKNESS,  SPACING  10.5  in.,  STEERING  20° 


In  each  of  Figs.  3  through  14,  the  ordinate  is  tne 
difference  frequency  pressure  with  dome,  normalized  to  the  dif¬ 
ference  frequency  pressure  without  dome.  The  abscissa  is  the 
dome  thickness  in  inches.  Each  figure  shows  the  results  obtain¬ 
ed  for  both  the  thick  and  thin  models .  In  Figs .  3  through  8 , 
the  result's  for  0°  steering  are  presented.  The  periodicity  of 

the  results  with  k  dcos0  is  very,  evident  in  these  results.  For 

P 

example.  Fig.  3  shows  good  agreement  between  thick  and  chin; 
furthermore,  this  separation  does  not  correspond  to  a  molt role 
wavelength  spacing  at  either  primary  frequenc>  in  the  dome- 
transducer  cavity;  therefore,  any  standing  wave  effects  are 
minimized.  On  the  other  hand,  in  Fig.  4  we  see  poorer  agree¬ 
ment  and  significantly  less  transmission  loss  than  in  Fig.  3. 

For  the  results  shown  in  Fig.  4  we  have  selected  a  dome- transducer 
spacing  resulting  in  a  resonance  condition  at  60  kHz.  This 
effect  is,  of  course,  amplified  by  the  ideal  geometry  and  lack 
of  damping  in  our  model. 

The  thick  model  resonance  out  around  1.0  inch  is  due 
to  a  quarter  wavelength  thickness  of  the  dome  at  one  of  the 
primary  frequencies.  Notice  that  this  internal  dome  resonance 
does  not  occur  in  Fig.  3  where  we  do  not  have  a  strong  standing 
wave  pattern  built  up  in  the  dome-transducer  cavity.  The  upturn 
of  the  thick  result  in  Fig.  3  is  caused  by  the  onset  of  a  half 
wavelength  resonance  in  the  material  which  occurs  at  a  thickness 
of  approximately  2  inches.  It  should  also  be  realized  that  the 
resonance  seen  in  Fig.  4  in  the  thick  model  is  very  sensitive 
to  dome -transducer  spacing  and  that  a  slightly  different  spacing 
would  reoult  in  the  disappearance  cf  the  peak  at  about  ]  inch 
and  merge  the  two  results  as  in  Fig.  3. 
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TABLE  II 


PARAMETERS,  MEASURED  DATA 


Primary  Frequencies 
(measured) 

Dome  Material 
h  (measured) 

p  (assumed) 

E  (assumed) 
v  (assumed) 
(assumed) 

Cg  (assumed) 
Fluid 

c  (assumed) 

p  (assumed) 

Dome- Transducer 
(measured) 

Steering  Angle 


246  kHz 
234  -kHz 

Steel 

0.087  in. 

.  lb --sec2 
7.35(I0)~4 

in. 

30(10)6  lb f/ in; 2 
0.3 

2.34(10)'>  in. /sec 
1. 25(10)5  in. /sec 
We  ter 

60,000  in. /sec 

-  lb --sec2 

9.66(10)"5  — - - 

in. 

3.0  in. 


Spacing 
(measured)  10.0 


1- 
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The  results  in  Figs.  5  through  S  are  like  those  in 
Figs.  3  and  4  but  for  different  dome- transducer  spaci:.gs.  The 
importance  of  these  result*,  is  that  they  show  the  periodicity 
of  ?-  with  half  wavelength  increases  in  dome- transducer  spacing 
Furthermore,  these  results  show  that  as  c  is  increased  the  dif¬ 
ference  in  cB  of  P„  between  a  resonant  and  antiresonant  dome- 
transducer  spacing  decrease. 

In  Figs.  9  through  14,  results  are  presented  for  a 

0 

steering  angle  of  20  .  For  these  off-axis  cases,  the  standing 
wave  patterns  are  different  than  for  d  =  0°  ;  hence,  we  do 
not  see  significant  effect. s  of  cavity  resonance,  for  the  same 
values  of  d,  as  we  did  in  Figs.  3  through  S.  The  main  effects 
of  changing  impedance  for  this  value  of  d  are  seen  for  small 
dome  thicknesses  and  are  prominent  in  Figs.  10  and  13. 

As  previously  stated,  we  were  aware  of  only  one 
piece  of  measured  data  available  for  comparison  with  our  model, 
and  this  comparison  is  shown  in  Fig.  15  with  appropriate 
parameters  given  in  Table  II.  Nothing  concrete  can  be  said 
about  our  models  or  their  predictive  capability  based  on  a 
comparison  with  one  measured  data  point;  however ,  it  was  encour 
aging  to  us  that  Konrad's  measured  point  was  bracketed  by  our 


two  models. 


CONCLUSIONS 


1.  Cavity  resonances  and  material  resonances  are 
important  effects  affecting  P_. 

2.  For  nonresonance  conditions  thick  and  thin  model 
results  agree  very  well  for  plate  thicknesses  much  greater  than 
A/10 . 
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3.  For  a  velocity  control  transducer,  if  the  imped¬ 
ance  the  transducer  sees  can  be  controlled  by  the  dome -transducer 
spacing,  then  one  can  achieve  more  power  out  with  a  dome  ..(under 

certain  conditions)  than  with  no-dome.  The  degree  to  which  this 

/• 

effect  would  carry  over  into  less  idealized  geometries  is  ques¬ 
tionable  however. 

4.  Our  two  models  bracketed  the  one  measured  data 
point  available  to  us. 

5.  The  following  calculating  procedure  is  suggested 
as  probably  being  adequate  to  make  predictions  of  the  influence 
of  a  dome  upon  the  performance  of  parametric  transmitting  array: 

1.  Utilize  existing  linear  dome  programs 

in  order  to  calculate  the  two  primary  fields 
at  the  dome's  outer  edge.  This  can  be  done 
both  for  thick  and  thin  dome  models. 

2.  Define  the  appropriate  Green's  function 
T  for  the  geometry  at  hand  which  satisfies 
the  boundary  condition  at  the  dome's  outer 
surface  of  T  -  0. 

3.  An  integral  similar  to  that  given  in 
Eq.  (25)  can  be  set  up.  This  will  lead  to 
the  definition  of  an  appropriate  coefficient 
A  which  will  be  sensitive  to  the  angle  of 
incidence  6  ,  the  frequencies  involved,  and 
the  dome  material  parameters.  Then  the  ratio 
of  A  (with  dome)  to  A  (without)  can  be  cal¬ 
culated  in  order  to  gain  insight  into  the 
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dome's  effect  upon  the  radiation  field  at  the 
difference  frequency.  It  should  be  emphasized 
that  for  three  dimensional  geometries  this 
approach  would  reveal  the  dome's  influence  upon 
side lobe  structure  in  both  the  primary  and 
secondary  pressure  fields  as  well  as  the 
efficiency  of  generation  of  P_ . 

7.0  RECOMMENDATIONS 

1.  Simple  experimental  measurements  need  to  be  made 
to  validate  this  model. 

2.  A  parameter  study  should  be; carried  out  using 
current  Navy  window  materials  to  identify  those  most  useful  for 
parametric  arrays. 

3.  The  concepts  of  this  model  should  be  analyti¬ 
cally  extended  to  more  realistic  geometries  and  experimentally  - 
verified  in  a  small  parallel  measurement  program. 
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